Work and heat probability distribution of an optically driven Brownian particle: 

Theory and experiments 



O 

o 

(N 

> 
O 

CO 



c3 



i 

C 

O 
o 



> 

On 

m 
q 

o 
l> 
o 



X 



A. Imparato, 1 L. Peliti, 2,3 ' 4 G. Pcscc, 2,5 G. Rusciano, 2,3 and A. Sasso 2 ' 3 

1 Dipartimento di Fisica and CNISM, INFN Sezione di Torino, 
Politecnico di Torino, C.so Duca degli Abruzzi 24, 10121 Torino, Italy 
2 Dipartimento di Scienze Fisiche, Universita "Federico II", 
Complesso Monte S. Angelo, 80126 Napoli, Italy 
3 CNISM, Napoli, Italy 
4 INFN, Sezione di Napoli, Italy 
5 CNR-INFM Coherentia, Napoli, Italy 
(Dated: September 15,2007) 

We analyze the equations governing the evolution of distributions of the work and the heat ex- 
changed with the environment by a manipulated stochastic system, by means of a compact and 
general derivation. We obtain explicit solutions for these equations for the case of a dragged Brow- 
nian particle in a harmonic potential. We successfully compare the resulting predictions with the 
outcomes of experiments, consisting in dragging a micron-sized colloidal particle through water with 
a laser trap. 
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The study of the physics of small systems has re- 
cently received a boost by the possibility of manipulating 
nanosystems and biomolecules. The fluctuations of the 
work and heat that these small systems exchange with 
the environment while being manipulated can be of the 
order or even larger than the thermal energy, leading 
to "transient" violations of the second principle of ther- 
modynamics. The distributions of heat and work have 
been experimentally studied for a few brownian systems 
§, |, |, |§. The probability distribution function (PDF) 
of the work done on a Brownian particle dragged by a 
moving quadratic potential was derived in Js], ^] . The dis- 
tribution turns out to be gaussian, what has been taken 
as an ansatz in |(| and confirmed in [Q by means of a 
rather involved path integral calculation. On the other 
hand, obtaining the PDF of the transferred heat repre- 
sents a much more difficult task: the Fourier transform 
of this function was obtained in refs. ^] by exploiting 
the energy balance and the gaussian ansatz for the work 
PDF, valid when the potential is quadratic. 

In the present paper we derive in a simple way the dif- 
ferential equations governing the evolution of the PDFs 
of the work and heat exchanged by a brownian particle, 
valid for any choice of the potential acting on the par- 
ticle. The solutions of these equation turn out to fulfill 
the well-known fluctuation relations. We evaluate the 
solution of these equations for a moving harmonic po- 
tential. We then experimentally study the work and the 
heat exchanged by a colloidal particle dragged through 
water by an optical trap. The PDF's predicted by our 
equations result in an excellent agreement with the ex- 
perimental data. We were inspired by the experiment of 
Wang et al. where the work done on a similar sys- 
tem was measured. However, in that experiment, only 
the performed work, and not the heat transferred, was 
sampled. Moreover, the expected gaussian distribution 
of the performed work was not verified, and a detailed 



comparison with the theoretical predictions was not at- 
tempted. However in a subsequent paper ||, the authors 
stressed that the PDF of the work has to be gaussian in 
their experimental conditions. 

Let us consider a Brownian particle in the overdamped 
regime, driven by a time-dependent potential U (x, X(t)), 
where X is an externally controlled parameter, that 
varies according to a fixed protocol X{t). The Langevin 
equation is given by 



UT 

^ = -ru'(x,x) + f(t), 



(i) 



where (/(*)/(*')) = (2l7/3)<5(t - t') and the prime de- 
notes derivative with respect to x. We have defined 
(3 = (fceT) -1 , and T = l/Qirnj, for a spherical particle 
with radius r, in a medium of viscosity rj. 

The thermodynamical work done on the particle is de- 
fined by 



(2) 



Besides work, the particle also exchanges with the envi- 
ronment heat, whose expression is 



= / dx^- = 

ox 



dt'x{t')U'(x(t'),X(t')). (3) 



If Q > the particle receives heat by the environment. 
Note that the integrals appearing in the last equation are 
stochastic integrals that must be interpreted according 
to the Stratonovich integration scheme ||. Let AU be 
the potential energy difference between the final and the 
initial state: the balance of energy for the manipulated 
particle reads AU = W + Q, which follows immediately 
from eqs. Note that the quantities Q and W have 

to be regarded as stochastic variables, whose value at 
time t depends on the specific stochastic trajectory. 
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The differential equation governing the time evolution 
of the PDF of the work is given by [pLOfj 



dt<l>(x,W,t) 



Td 2 (j) 
^dx 2 



X 



dU 



dX dW 



(4) 



It can be easily shown that the solution of eq. (GJ) satisfies 
the Jarzynski equality |lC| ]. 

The differential equation for the joint PDF ip(x, Q, t) of 
the position x and the heat Q is obtained as follows. In a 
short time interval 8t the heat exchanged by the particle 
with the environment reads SQ = dU — dtU St = U'dx, 
and thus the time derivative of Q is given by 

^=U'^=U'(-TU' + f(t)). (5) 

Thus equations (Q) and (||) describe two coupled stochas- 
tic processes. We now define the vectors of the stochastic 
variables y and of the forces F by 



and the diffusion matrix 



-TU 1 
-TU' 2 



B 



r/& (T/p)u' 



- - \ (T/(3)U> {T/p)U r - 



(6) 



(7) 



Then the differential equation governing cp(x,Q,t) 
straightforwardly follows [JO} [l2| : 

M*.0,i)=~(lV) + £(B.A^. (8) 

By introducing the generating function x(ir,A,i) = 
J dQ exp(AQ)(p(x, Q, t), we obtain the simpler equation 

that was first derived by Lebowitz e Spohn [[u|. Note 
that the operator appearing on the rhs of this equation 
changes into its adjoint by the substitution A — ► (3 — 
A. As discussed in ref. jl]|, this symmetry implies the 
Gallavotti-Cohen fluctuation relation (l4| for our system. 

By defining the function g(x, A, t) as x{x, A, t) = 
g(x,X,t) exp [-6(X)U(x,t)/2], with S(X) = /3-2X, eq. @ 
becomes 



d t g 



(3dx 2 



U' 2 
TP—g- 



\u" 9+ 5 -^gd t U. (10) 



This equation has the form of an imaginary-time 
Schrodinger equation. 

It is worth remarking that eqs. (EJ), and (||) hold 
for any choice of the potential U(x, X{t)). Moreover, the 
present approach can be easily generalized to the case of 
a Brownian particle with inertia H, 15, fl6| . 



We now consider the particular case of a harmonic po- 
tential 



U{x,t) = -{x-X{t))\ 



(11) 



with X{t) = vt, i.e., the center of the potential moves 
with a constant velocity v. We shall assume that the 
particle is initially in thermal equilibrium, with the po- 
tential centered at X — at t — 0. It is then possible to 
solve directly equation A), obtaining 



(x, W, t) = Af t exp 



W-W(x,t)) g k 

^(t) -T <*-««)>' 

(12)" 

In this equation, having defined t = 1/Fk and a(t) = 
e~*/ T , we have £(*) = vt (a(t) - 1 + t/r), a 2 {t) = 
v 2 T 2 kf3- 1 [2t/r + 1 - (2 - a(t)) 2 }, 

W(x,t) = tv 2 rk(2- a(t)) - vxrk(l - a(t)) (13) 
-v 2 T 2 k(2 + a 2 (t)-3a(t)); 

A/T 1 = ^(ttvt/P) 2 (2t/r + 1 - (2 - a(i)) 2 )(.14) 

The unconstrained PDF $(W, t) = f dx <p{x,W,t), is 
given by 



$(W,i) = A"/ exp 



(W-u 2 r 2 fc(a(i)-l + f/r)) 2 
4v 2 r 2 fc(a(f) - l+t/r) 



(15) 

where Af( = [An/3' 1 v 2 T 2 k (a{t) - 1 + t/r)] 1/2 . A simi- 
lar result was first obtained in g in a special case, and 
then in ||, by using qualitative arguments and by as- 
suming that <&(W,t) is gaussian, and more recently in 
@' by a functional integral technique. We now see that 
eqs. (|l2 ,P~5[) can be straightforwardly derived as solutions 



of eq. (4). An approach analogous to ours was used in 
ref. fll7|l7 leading again to eq. (^5|). 

We now turn to the heat PDF: substituting eq. ( |Tl| ) 
into eq. ( |To| ) one obtains a Schrodinger-like equation for 
the harmonic oscillator, that can be solved exactly. As- 
suming that the particle is at thermal equilibrium at 
t = 0, with v = for t < 0, the solution of eq. (||) 
reads 



x(x,X,t) = exp 



.& U(X)X (t))-^z(x,t) 



^ e lnt c n (X)il; n (z(x,t)), 



x y e 

n=0 



(16) 



where j n = (-n/r + 5 2 {\)v 2 /AT(3 - (3v 2 /AT), and 
z(x,t) = x — vt + 5(X)vt/P, and where ifj n (z) are the 
eigenf unctions of the Schrodinger equation for the har- 
monic oscillator, with the substitutions h 2 /m — > 2T//3 
and mux 2 — > Tfjk 2 /2. The value of coefficients c„(A) is 
determined by the initial condition x{ x ,X,t — 0). Note 
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that r = 1/rfc sets up the characteristic time scale for 
both work and heat fluctuations. 

Case a): v — for t > 0, i.e., a fixed potential. The be- 
havior of the generating function \&(A, t) = J c&c x(^7 A, i) 
in the long-time limit is governed by the eigenfunc- 
tion tpo(z(x,t)) associated with the smallest eigenvalue. 
Thus, after some algebra, one finds ^S(X,t — > oo) = 
1/^/1 — (A//3) 2 . Therefore, in the case of constant po- 
tential, the heat unconstrained PDF in the long time 
limit has the expression 



-»-/ 



dX 



T *(A,t^oo)e" AQ = /3 



MP\Q\) 



<p(Q,t- - , • 

J Iiri vi 

(17) 

where Kq{x) is the zero-th order modified Bessel function 
of the second kind. 

Case b): v > 0. Also in this case the long time be- 
havior of the solution of eq. (||) will be dominated by the 
eigenfunction with n = 0, which is a gaussian function. 
Thus, in the long time limit, one finds 



exp 



4r 



+ 



-> oo) 
2A(3-4(A//3) 2 ) 



4A 

J 



Tk(f3 + \) 



-I t 



(18) 



-1/2 



and, integrating by the saddle-point method, one finds 



(p(Q, t — > oo) 



1 

2?ri 



dA*(A,i-> oo)e AC? (19) 



exp 



r/3 



(20) 



Note that in the long time limit IF + Q = 0. As discussed 
m ref. §, some care has to be taken when calculating the 
integral fll9| ) with the saddle point method. The resulting 
calculation shows that the function (piQ, t) is gaussian up 
to a subleading term of order \j\ft 0. 

In order to test these results, we have experimentally 
observed the trajectories of a colloidal particle in an opti- 
cal trap, which is well described by a quadratic potential 
( pi] ) near its focus X(t). 

The Optical Tweezers system consisted of a home 
made optical microscope with a high numerical aper- 
ture water immersion objective lens (Olympus, UP- 
LAPO60XW3, NA=1.2) and a frequency and amplitude 
stabilized Nd-YAG laser (A = 1.064 fj,m, 500 mW, In- 
nolight Mephisto). The sample cell was made with a 
glass coverslip of 150 /im thickness and a microscope slide 
glued together by parafilm stripes of about 100 fim thick- 
ness. Polystyrene micro-spheres produced by Postnova 
(density: 1.06 g/cm 3 , refractive index: 1.65) with a di- 
ameter of 2.00±0.05 /xm were diluted in distilled water 
to a final concentration of about l-=-2 particles/pl. The 
sample cell was mounted on a closed-loop piezoelectric 
transducer stage (Physik Instrumente PI-517.3CL) which 
allowed movements with nanometer resolution. Moving 
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FIG. 1: Histogram of the work exerted on the colloidal parti- 
cle by the optical trap, as given by eq. (Q, for t = 0.01s<r 
(main figure) and f = 0.5s > t (inset). The lines correspond 
to the expected function (|l5|), with no adjustable parameter. 
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FIG. 2: Histogram of the heat exchanged by the colloidal 
particle with the environment, with fixed optical trap v = 0, 
as given by eq. ([]), for t = 0.5 s. The line correspond to 
the expected PDF as given by eq. (|l7|), with no adjustable 
parameter. Inset: same data, with logarithmic ?/-axis. 



the stage in a given direction corresponds to moving the 
optical trap focus in the opposite direction. The sample 
temperature was not stabilized but continuously moni- 
tored using a negative temperature coefficient thermis- 
tor positioned on the top surface of the microscope slide. 
The temperature during a complete set of measurements 
remains constant T = 296.5 K, within 0.2 K. The trapped 
bead was positioned in the middle of the sample cell to 
avoid any surface effects. The thermally driven motion 
of trapped beads was monitored by a InGaAs quadrant 
photodiode (Hamamatsu G6849) placed in the back fo- 
cal plane of the condenser lens [Q. The response of 
our quadrant photodiode was linear for displacements of 
about 300 nm with a resolution of 2 nm, and its band- 
width was about 250 kHz. The trajectories in the trans- 
verse x-y plane were sampled at 125 kHz using a digital 
oscilloscope (details on the experimental setup can be 
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FIG. 3: Histogram of the heat exchanged by the colloidal par- 
ticle with the environment, with v = 1 fim/s, and t = 0.5 s 3> 
r, as given by eq. (^|). The line correspond to the expected 
PDF as given by eq. (M), with no adjustable parameter. 



found elsewhere Jl9|]). The duration of each trajectory 
measurement was 10 s. During the first 5 seconds the 
stage was at rest and we use this period to compute the 
power spectral density of the particle position to obtain 
the trap stiffness, which takes the value k — 6.67 x 10~ 7 
N/m, and the calibration factor of the quadrant photo- 
diode 1 20 . Then at t — 5 s the stage started to move 
with a speed of v = 1/im/s for 5 s along one axis. After 



a pause of 1 second the above described sequence started 
again, but the stage was moved in the opposite direction. 
We repeated this procedure 300 times and back and forth 
trajectories were recorded for further analysis. Thus the 
overall number of trajectories considered is 600. Note 
that for each trajectory we measure both the work done 
on the particle and the heat, as defined by eq. (0) and 
(||), respectively. We have, under our experimental condi- 
tions, T = 1/67T7-77 = 5.76 x 10 7 s 2 /kg, yielding r ~ 0.026 
s. 

In figure |l|, we plot the histogram of the work ex- 



erted on the particle by the optical trap, finding a good 
agreement with the expected PDF Q(W, t), as given by 
eq. (|l5|): the distribution of the work turns out to be 
gaussian both at short and long times. At short times 
(t = 0.01 s) the gaussian is peaked around zero, but the 
peak moves to positive values of W as t increases. 

In order to evaluate the heat Q from the particle tra- 
jectories, we exploit the discrete version of eq. (||): Q — 
E t< K - xti-r) k/2 [{x ti - X{U)) + {x u _ x - Xik-i))] . 
We plot the histogram of the measured heat for the 
motionless trap, with t — 0.5 s, in fig. |[ The histogram 
agrees nicely with eq. ( |l7| ) , in particular the tails of the 
distribution are found to be exponential. This behavior 
was also found at shorter times (data not shown). 

Finally, in fig. (||) the histogram of the measured heat, 
for the trap moving with v = 1 /im/ s, is plotted, in the 
long-time range t — 0.5 s»t. The distribution is found 
to be gaussian, in agreement with eq. (|2p|). By comparing 
figures and |, it can also be seen that the mean values 
of the work and of the heat are the negative of each other, 
as expected in the long-time range. At shorter times we 
observe that the tails of the distribution of the measured 
heat fall off exponentially, with time-dependent slopes 
(data not shown). We have noticed that for observation 
times longer than 0.5 s both the heat and work distri- 
butions appear slightly broader than the theoretical pre- 
dictions, whereas their centers remain in good agreement 
with the expected ones (data not shown). We ascribe 
this fact to the presence of low-frequency (smaller than 
1 Hz) noise affecting our experimental set-up. 

We have shown that it is possible to solve explicitly the 
differential equations governing the evolution of the PDF 
for the work and heat exchanged by a dragged brownian 
particle, and that the resulting predictions are vindicated 
by experiment. In particular one observes a non negligi- 
ble probability for the "transient violations" of the sec- 
ond law of thermodynamics, i.e., positive values of the 
exchanged heat Q. 
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